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SKIN  TEMPERATURE  RISE  IN  AN  AIRCRAFT 
EXPOSED  TO  THERMAL  RADIATION 
FROM  A  NUCLEAR  EXPLOSION 


SUMMARY 


It  is  shown  that,  in  mnny  siyoations,,  aei'odynajiic  cooling  substantially 
reduces  the  tempera t^are  rise  in  thin  skins  of  edrcraft  exTx>sed  in  flight  to 
thermal  radiation  from  nv\clear  explosions.  The  maxiinum  acceptable  thermal 
dose,  where  this  is  determLied  by  skin  teii5>5rature  rise,  depends  upon  the 
time  scale  of  the  theriLal  input  aixi  thus  upon  the  yield  oi*  the  weapon 
involved. 

The  problem  is  attacked  by  a  finite  difference  analysis,  and  the 
maximum  tcu^ierature  rise  and  the  time  at  v-hich  it  occurs  are  sham  to  depend 
on  a  single  non-d.iniensional  jjarameter.  The  results  are  presented  for  a 
range  of  this  parameter  corresponding  to  the  majority  of  practiool  situations. 


1  INTRCDUCTION 


In  many  situations,  the  maximum  yield  of  weapon  which  an  aircraft  may 
safely  deliver  is  governed  by  the  thermal,  dose  it  can  withstand.  In  parti¬ 
cular,  the  temperature  rise  in  the  skin  may  be  limited  to  some  critical  value 
determined  by  deterioration  of  the  properties  of  the  material  or  by  unaccept¬ 
ably  high  thermal  stresses.  It  is  clear  that  the  operational  capabilities 

of  an  aircraft  can  only  be  detemined  realistically  if  both  the  critical 
temperature  limits  and  the  temperature  rise  in  a  given  situation  are  known 
accurately.  This  Note  considers  the  latter  factor. 

In  previous  British  analyses  of  the  effects  of  theriiial  radiation  on 
aircraft,  it  has  been  assLined  that  the  incident  thermal  dose,  corrected  for 
reflection  at  the  sui'face  and  the  obliquity  of  the  skin  to  the  direction  of 
propagation,  is  thereafter  totally  absorbed  by  the  skin  and  contributea  in 
its  entirety  to  the  temperature  rise.  It  \/as  known  that  during  the  tinK  in 
which  the  thermal  zad.iation  was  received,  soTiie  heat  would  be  lost  from  the 

surface  by  forced  convection  to  the  airstream,  but  this  effect  had  not  been 

evaluated  and  had  been  considered  to  be  a  factor  of  safety.  With  the  advent 
of  thermonuclear  weapons  of  high  yield,  it  has  now  become  necessary  to 
examine  this  and  other  safety  factors,  in  order  to  determine  whether  these 
weapons  can  be  delivered  at  all  by  existing  aarcraft  without  damaging  them 
seriously. 
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This  Note  considers  the  effect  of  aerodyn4>/nic  cooling  on  the  temperature 
rise  of  aircraft  skins  exposed  to  thermal  radiation  from  nuclear  weapons. 

The  analysis  is  generalised  end  the  range  of  variables  chosen  is  believed  to 
cover  any  likely  situation# 

2  RATE  OP  THERMAL  INPUT 


The  thermal  energy  from  a  nuclear  explosion  is  released  in  a  pulse, 

the  rate  of  output  rising  rapidly  to  a  peak  and  decaying  more  slowly  over  a 

relatively  long  period  of  time.  It  is  possible  to  represent  the  thermal 

pulse  in  a  generalised  form  giving  the  rate  of  thermal  energy/  releeise  as  a 

fraction  of  its  peak  valvie,  in  terms  of  a  non-dimensional  time  whose  units 

are  the  time  t^^^  at  which  the  energy  release  is  a  maximum,  i.e. 
max 

p(t)/p  =  f(t/t  )  .  (i) 

^  **  max  ^  '  max'  ' 

1 

This  function,  which  is  independent  of  yield,  is  shov.n  in  Pig.1  ;  it  also 
represents  the  chape  cf  the  pulse  received  by  a  surface  exposed  to  the 
radiation. 

In  the  subsequent  analysis,  the  emoxmt  of  energy  received  up  to  a 
given  time  is  requixed,  aM  this  is  obtained  from  an  integral  of  equation  (l) 
and  takes  the  non-diiusnsional  form 

Q(t)/Q  =  f(t/t^)  (2) 

where  Q  is  the  total,  thermal  energy  received.  A  curve  of  equation  (2),  of 
American  origin,  has  been  widely  published^  end  is  shown  in  Pig. 2.  The  data 
actually  ;ised  in  thi.s  study  are  very  sirilar  to  those  represented  in  Pig.  2, 
but  were  obtained  in  British  trials  with  w'eapons  in  the  megaton  range^# 

There  appear  to  be  no  data  for  times  longer  than  10  t-^^y,  and  it  is  of 

importance  to  note  that  only  abo*at  of  the  theimal  energy  is  released  up 

to  this  time.  Since  j_ 

t  :=  0.032  secs  (Ref.  1)  (3) 

max  '  '  '  ' 


it  may  be  seen  that  some  20^5  of  the  theiroal  energy  is  still  to  be  released 

after  a  time  t  s  0. 32  W^,  which  for  a  5  Mt  weapon,  for  example,  is 
22.6  secs.  The  fact  that  most  of  the  energy  is  released  after  the  peak 
peimits  a  significant  amount  of  convective  cooling  to  occur  in  most  circum- 
stanoes,  as  shown  below. 

3  RATE  OP  HEAT  LOSS 


It  is  assumed  that  the  rate  of  loss  of  heat  from  an  aircraft  skin  to 
the  oirstream  can  be  represented  by  a  heat  transfer  coefficient  h  such 
that  ^en  the  temperature  difference  between  the  skin  and  the  airstream 
is  T,  the  rate  of  heat  loss  per  unit  area  is  hT.  It  is  not  prox^osed  to 
enter  into  a  detailed  examination  of  means  by  which  the  heat  transfer 
coefficient  may  be  calculated  in  different  circumstaiices,  but  a  brief  state¬ 
ment  of  relevant  data  is  given  below. 

The  heat  transfer  coefficient  at  any  point  depends  upon  the  conditions 
in  the  boundary  layer  and  the  location  of  the  point.  It  is  usually 
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assumed  that  the  development  of  the  boundary  layer  is  similar  to  that  on  a 
flat  plate,  for  which  considerable  experimental  data  are  available#  A 
number  of  empirical  formulae  have  been  developed  to  represent  these  data, 
and  those  recommended  for  use  in  kinetic  heating  studies  seem  to  be  the 
most  appropriate^.  At  any  point  whose  distance  is  x  downstream  from  the 
leading  edge,  the  local  heat  transfer  coefficient  is  given  by 

h  =  0,332  pVCp  Pr  Re  ^  (4) 


when  the  flow  in  the  boundaiy  layer  is  laminar,  and  by 

h  t,  0.176  pVCp  (log^Q  (5) 

when  the  flow  in  the  boundary  layer  is  turbulent#  These  equations  may  be 
modified  by  compressibility  effects  and  chordwise  pressure  gradients  and 
the  heat  transfer  coefficient  also  varies  somewhat  with  the  temperature 
difference  T  when  this  becomes  large#  The  effects  of  soiae  of  these  factors 
are  oonsidered  elsewhere,  and  appropriate  specialist  reports  should  be 
consulted3>4,5.  It  will  be  seen  later,  however,  that  the  maxijiium  tempera¬ 
ture  rise  in  many  aituations  is  relatively  insuisitive  to  the  heat  transfer 
coefficient,  to  an  extent  which  makes  great  accuracy  in  determining  the 
latter  uxmecessary# 

4  RATE  OP  TEMEERATURS  RISE 


In  attempting  a  general  analysis  of  the  rate  of  increase  in  skin 
temperature,  it  is  necessary  to  meto  the  following  assumptions,  all  of 
which  appear  to  be  justified  in  the  majority  of  situations: 

(a)  The  heat  transfer  coefficient  reanoains  constant# 

(b)  The  themal  properties  of  the  skin  material  remain  constant. 

(c)  The  reflectivity  of  the  surface  remains  constant# 

(d)  The  skin  and  boundary  layer  ten^jeratures  are  initially  equal# 

(e)  There  are  no  temperatiure  gradients  through  the  skin. 

(f)  Heat  loss  to  adjacent  structure  is  negligible# 

(g)  Heat  loss  by  radiation  is  negligible. 

Under  these  circumstances  the  ten^^erature  rise  in  the  skin  is  governed 
by  the  differential  equation 


zcl  ^  =  P(t)  -  hi 


(6) 


where,  in  this  context,  P(t)  is  the  rate  of  heat  input  per  unit  area  of  skin 
at  time  t. 


Cleai'ly,  equation  (6) 
analytic  function  of  time. 


con  only  be  evaluated  if  P(t)  is  known  as  an 
Since  it  is  not,  it  is  necessary  to  consider  an 
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ane lysis  using  finite  differences  of  the  variables*  Equation  (6)  then 
becomes 


or 


AT  = 


■  S  =  p(t)-hi 

(7) 

P(t)  At  hT  At 
""  zcA  zcA 

(8) 

where  T  is  taken  as  the  temperature  rise  at  the  centre  of  t^e  time  interval 
At.  ^ow  P(t)  At  is  the  amount  of  heat  entering  the  skin  during  the  interval 
At,  which  is  AQ(t).  Then 


'  zct  zat  • 


(9) 


To  effect  a  general  solution,  we  may  make  tlie  substitutions 


and 


Further,  let 


dQ(t)  ,  Q 


At 


At 

t 

max 


^  • 

max 


2cA 


T 


i 


(10) 


where  T.  is  the  "ideal"  ten^jerature  rise  which  would  resxilt  if  no  heat  v/ere 
lost,  and  let 


ht 


max 

zcA 


b 


(11) 


where  b  is  a  non-dimensional  variable* 
Then  equation  (9)  becomes 


whence 


AT 


At 

t 

max 


T 


T 

^i 


Q 


At  2. 

t  T. 
max  1 


a  non-dimensional  equation* 


(12) 


(13) 
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Then  during  the  nth  interval  of  time  the  non-dimensional  temperature 
rise  increases  by 


AQ(t)^ 

Q 


At 

•  t 

max 


T 


n 


(14) 


At  the  end  of  this  interval,  let  the  non-dimensional  ten5)erat\ire  rise 
be  —  ,  where 


IT"  T  AT 

n  _  _n  1  _ n 

rp  “  rp  •  rp 

^i  ^i  1 


(15) 


During  the  next  interval,  the  non-dimensional  ten^jerature  rise  increases 
by 


But 


Q 


b 


(16) 


AT 


n+1 


T. 

1 


(17) 


Then  equation  (16)  gives 


Now  the  non-dimensionaJ.  tenqperature  rise  at  the  end  of  the  (n+1)th 
interval  is  given  by 


AT 


T. 

1 


(19) 


and  substituting  for 


from  equation  (17) 


(20) 
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It  is  now  possible  to  erect  a  series  from  which  the  teinperature  rise  at 
any  time  can  be  determined.  During  the  first  iriterval 


(21) 


and  thus,  f rom  e  quation  (llf) 


and 


(22) 


(23) 


Subsequently 


and 


T 


i 


(20) 


Using  equations  (22),  (23),  (18)  and  (20)  in  turn,  the  variation  of 
T/T^  with  determined  for  any  value  of  the  non-dimensional 

variable  b»  Thus  the  ten^rature  rise  time-history  for  any  situation  may  be 
presented  in  terms  of  cue  parameter.  In  using  this  method,  it  should  be  noted 
that  the  temperature  rise  occurs  at  the  centre  of  the  nth  interval,  i.e.  at 

a  time  (n— 5)At*  In  order  to  produce  a  generalised  solution,  a  large  number 

of  temperature -time  series  has  been  evaluated,  using  the  basic  assumptions 
given  previously.  It  is  worth  noting  however,  that  in  any  particular 
situation,  the  solutions  may  be  obtained  if  all  the  assun5)tions  except  (e) 
are  discarded,  so  long  as  the  variation  of  all  the  relevant  terras  is  known  as 
a  function  of  temperature  or  of  time. 

Typical  solutions  have  been  worked  out  in  the  range  0.01  <  b  <  1.0, 

within  which  most  practical  situations  are  believed  to  lie.  Some  of  these 

are  shown  in  Pig.  3,  the  working  having  been  condiided  just  after  the  mrjcima 

had  been  reached.  For  these  solutions,  the  time  intervals  used  were 

At/t  =  0. 25. 

'  max 

In  Pig3«/f  and  5,  the  inaxiicum  non-dinienaional  ten^erature  rise  T 
and  the  non-dimensional  time  which  it  occiirs  are  plotted  against 
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the  parajneter  bn  A  difficulty  was  encountered  in  determining  these  values 
when  they  occurred  after  10  since  there  are  no  input  data  beyond  this/'/ 

point,  y&iich  corresponds  approximately  to  b  =  0.01  •  In  Fig«4>  the  extra¬ 
polation  down  to  b  =  0«005  has  been  dotted  to  eii5)haLSise  its  hypothetical 
nature.  This  extrapolation  is  made  since  situations  involving  low-yield 
weapons  and  little  cooling  have  been  found  to  involve  values  of  b  some¬ 
what  below  0*01  • 

5  DISCUSSION 

It  may  be  seen  from  Fig#4  that  the  maximum  ten^rature  rise  in  an 
aircraft  skin  in  the  presence  of  aerodynamic  cooling  is  substantially  less 
than  has  been  assumed  hitherto  when  no  allcAvance  for  cooling  has  been  made. 
The  variation  of  maximum  ten^rature  rise  with  the  parameter  b  is  not  rapid 
over  the  range  0.01  <  b  <  0.1,  which  is  believed  to  bo  the  most  important 
range.  It  therefore  appears  that  the  values  cf  the  properties  ^ich 
compose  b  need  not  be  known  v/ith  extreme  accuracy  in  order  to  determine  the 
maximum  temperature  rise  with  some  confidence. 

In  using  Fig.4,  it  must  be  remembered  that  although  high  values  of  h 
and  t^^^^  invariably  favour  low  temperature  rises,  small  values  oi  (zo^)  not 

only  increase  the  value  of  b,  but  also  that  of  T^,  since 


(10) 


Thus,  thin  skins  always  suffer  greater  temperature  rises  than  thicker  ones, 


This  is  illustrated  in 


even  though  they  favour 


Fig.6,  for  example,  where  the  actual  maximum  temperature  rise  is  plotted 

against  ht  for  skins  of  various  thicknesses  in  DTD. 5^6  material.  Fig.6 
^  max  ® 

also  demonstrates  more  clearly  than  does  Fig. If  the  relative  insensitivity 

of  the  maximum  temperature  rise  to  the  heat  transfer  coefficient.  Even  in 

the  thinnest  material,  a  change  in  h  by  a  factor  cf  two  changes  T  by  about 

10-2095. 

When  the  skin  is  very  thin,  the  reflective  paint  applied  to  the 
surfaces  may  contribute  significantly  to  the  heat  capacity.  If  the  seme 
assumptions  are  rnede  as  in  Section  subsequent  analysis  is  applicable 

when  the  product  (zc^)  is  replaced  by  the  total  heat  capacity  per  unit  area 
of  the  akin  together  with  its  paint. 

It  is  important  to  note  that,  when  aerodynamic  cooling  is  taken  into 
account,  it  is  no  longer  possible  to  state  a  maximum  tliermal  dose  which 
a  given  aircraft  will  withstand,  when  this  is  limited  by  the  rise  in  skin 
temperature.  Any  .such  critical  dose  must  be  associated  with  the  yield 
of  the  weapon  from  which  it  emanates*  In  general,  a  greater  total  themal 
dose  may  be  accepted  from  a  weapon  with  a  high  yield  than  from  one  with  a 
low  yield,  because  cf  the  increase  in  the  time  scale  for  the  former  weapon. 


6  ILTUSTRATIVE  EXAMPLE 

An  example  involving  a  hypothetical  situation  is  worked  out  below  as 
an  illustration  of  the  use  of  the  results. 

The  maximum  temperature  rise  in  the  aileron  skin  of  a  bomber  is 
required  when  a  total  thermal  dose,  corrected  for  reflectivity  and  obliquity 


11 

CONFIDENTIAL 


of  tne  surface,  of  10  cals  cnT^  is  received  from  a  2  Utm  explosion.  The 
skin  is  26  gauge  light  alloy  to  DTD  5^-6  and  the  point  of  interest  is  12  ft 
downstream  from  the  leading  edge.  The  aircraft  is  flying  straight  and 
level  at  M  =  0.8  at  45,000  ft. 

It  will  be  assumed,  in  this  example,  that  the  temperature  rise  is 
small  enough  for  the  properties  of  the  air  to  be  evaluated  at  ambient 
conditions,  i.e., 

p  =  4.6  X  10*^  slugs  X  ft"*^ 

^  -1  -1 
|i  s  2.54  X  10  slugs  X  ft  sec 

X  B  12  ft 

V  =  795  ft  seo"'’ 

Then  Re  b  ®  1*49  x  10^.  The  homdary  layer  would  thus  he  t\irhulent, 

and  froQ  equation  (5). 

h  s  0.176  X  pVCp  (log^Q  Re)"^*'^^ 

h  s  3.29  X  10“^  CHU  ft"^  aec""'  oq"'’ 

=  1.61  X  10**^  cals  cm”^  sec  ^  ®C  ^ 

Prom  equation  (3) 

1 

t  =  0.032  seo3 

max 

=  1.43  secs. 

For  26  gauge  nrD.546 


(zcC)  =  0.0262  cals  cm"^  oc"”'. 
Then  from  equation  (ll) 


h 


ht 


max 

zcl 


=  0.088 


Prom  Pig.4 


T^Ti  =  0.56 
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Prom  equation  (10) 


Then 


Prom  Pig.  5 


Then 


=  Q/zot 

B  381°C 

T  =  213°C  . 
in 


max 


3.75  . 


t  =  5«36  secs, 

m 

Thus,  the  maximum  t emperature  rise  is  213^0  occurring  5o36  secs  after 
burst. 


7  CONCLUSIONS 


A  simplified  analysis  has  been  made  of  the  effects  of  aerodynamic 
cooling  on  the  temperature  rise  of  thin  skins  exposed  to  thermal  radiation 
from  nuclear  weapons. 

It  is  shewn  that  the  maximum  temperature  rise,  expressed  as  a  fraction 
of  the  value  which  would  result  v/ithout  cooling,  is  a  function  of  a  single 
non-dimensional  parameter.  A  curve  is  given  of  this  function  together 
with  a  curve  shewing  the  time  at  which  the  maximum  ten^srature  rise  occurs, 
also  in  non-dimensional  fonn# 

The  results  indicate  that  the  maximum  ten^rature  rise  is  substan¬ 
tially  reduced  by  aerodynamic  cooling  in  many  practical  situations. 

Attention  is  drawn  to  the  fact  that  the  maximum  thermal  dose  which  a 
structure  may  withstand  depends  upon  the  yield  of  the  weapon  frexe  which 
it  emanates. 


NQTjftTION 

The  symbols  used  are  defined  as  follows: 
ht 

.L  ^  max 

b  parameter 

c  specific  heat  of  skin  material 

Cp  specific  heat  of  air  at  constant  pressure 

h  heat  transfer  coefficient 

k  thoimal  conductivity  of  air 
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NOTATION  (CONTD.) 


t 

skin  thickness 

t 

tiioe  after  burst 

t 

max 

time  at  which  P(t)  =  P 

'  '  max 

t 

m 

time  at  which  T  =  T 

m 

At 

interval  of  time 

X 

distance  of  point  from  leading  edge 

z 

density  of  skin  material 

Pr 

Prandtl  nunber 

P(t) 

rate  of  emission  or  reception  of  themal  energy  at  time  t 

^max 

maximum  value  of  P(t) 

Q(t) 

thermal  energy  emitted  or  received  up  to  time  t 

AQ(t) 

increment  of  Q(t)  in  time  interval  At 

Q 

total  thermal  energy  emitted  or  received 

Re 

X 

Reynolds  nunibe  r 

T 

temperature  rise 

"ideal**  temperature  rise 

T 

m 

maxirraom  temperature  rise 
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temperature  rise  during  nth  interval 
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weapon  yield  in  kilotons 

P 

density  of  air 

P 

viscosity  of  air 
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